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ABSTRACT 

Just as correlations between fluctuating radial and azimuthal velocities pro- 
duce a coherent stress contributing to the angular momentum transport in tur- 
bulent accretion disks, correlations in the velocity and temperature fluctuations 
produce a coherent energy flux. This nonadvective thermal energy flux is always 
of secondary importance in thin radiative disks, but cannot be neglected in non- 
radiative flows, in which it completes the mean field description of turbulence. It 
is, nevertheless, generally ignored in accretion flow theory, with the exception of 
models explicitly driven by thermal convection, where it is modeled phenomeno- 
logically. This flux embodies both turbulent thermal convection as well as wave 
transport, and its presence is essential for a proper formulation of energy conser- 
vation, whether convection is present or not. The sign of the thermal flux is likely 
to be outward in real systems, but the restrictive assumptions used in numerical 
simulations may lead to inward thermal transport, in which case qualitatively 
new effects may be exhibited. We find, for example, that a static solution would 
require inward, not outward, thermal transport. Even if it were present, thermal 
convection would be unlikely to stifle accretion, but would simply add to the 
outward rotational energy flux that must already be present. 

Subject headings: accretion — accretion disks — black hole physics — instabilities 
— (magnetohydrodynamics:) MHD — turbulence 



1. Introduction 

The fate of gas in the vicinity of a black hole remains a topic of intense astrophysical interest. 
The classical motivation for black hole accretion is of course the possibility of prodigious 
energy output (Lynden-Bell 1969), but in more recent years attention has also been drawn 
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to a class of flows that radiate inefficiently (Abramowicz, et al. f988; Narayan, & Yi f995; 
Narayan, Igumenshchev, & Abramowicz 2000). With the unambiguous Chandra finding 
that the Galactic Center source Sgr A* is severely underluminous (Baganoff, et al. 200f), 
these models have now acquired a compelling observational motivation. 

Analytic models of nonradiating accretion flows (hereafter NRAFs) are generally based 
on the classical a prescription, whereby all transport is subsumed into a single enhanced vis- 
cosity parameter. NRAFs are very likely to be turbulent, of course, because the combination 
of differential rotation and a magnetic field is prone to the magnetorotational instability, or 
MRI (Balbus & Hawley 1991). The simplifying assumption that has almost universally been 
made is that the sole dynamical effect of this turbulence is to act as the enlarged a viscosity 
parameter. In earlier NRAF models, this assumption applied both to energy as well as to 
angular momentum transport (Narayan, & Yi 1994; Blandford & Begelman 1999). More 
recently, a class of accretion flow has appeared that specifically invokes thermal convection 
as the dominant mode of energy transport, while at the same time arguing that the angu- 
lar momentum stress tensor effectively vanishes (Narayan, Igumenshchev, & Abramowicz 
2000). 

The purpose of this paper is to elucidate, from first principles, the nature of turbulent 
transport couplings in an NRAF. The formalism that allows a connection to be made be- 
tween phenomenological transport modeling and the fundamental fluid equations is weak 
turbulence theory. This is a less restrictive limitation than it may sound. It is equivalent to 
retaining the quadratic correlations in the turbulent fluctuations, and is a mathematically 
self-consistent approach. Phenomenological models of turbulent accretion flows certainly do 
not attempt higher order accuracy, and often settle for far less. 

Applications of this type of formalism to nonradiative flows were carried out by Quataert 
& Narayan (1999), and to black hole accretion more generally by Kato, Fukue, & Mineshige 
(1998). But in its emphasis upon modes of transport rather than difficult questions of closure, 
our approach more closely follows that of Balbus & Papaloizou (1999). These authors used 
weak turbulence theory to elucidate the foundations of a modeling in thin radiative disks, 
and found that whereas MHD processes naturally lend themselves to such a description, 
self-gravity generally does not. A related finding of this study was that the only correlation 
tensor of dynamical importance that emerges in an MHD turbulent Keplerian disk is the 
Maxwell- Reynolds stress. This fact, in turn, can be traced to the assumption that the ratio of 
the isothermal sound speed cs to Keplerian rotation velocity vk satisfies cs/vk <C 1- When 
this assumption no longer holds, other transport coefficients become dynamically significant 
(as happens in self-gravitating disks), in which case the simple viscous model itself breaks 
down. 

The regime Cs/vk ~ 1 coincides with the NRAF regime. We are led, therefore, to 
consider flows in which the energy extracted from the differential rotation is, in addition 
to being locally dissipated, nonlocally transported by other correlations associated with a 
turbulent or wavelike energy flux. This form of transport has the potential to be important 
in any NRAF, whether or not thermal convection is associated with the accretion flow. Its 
role is quite distinct from that of the dynamical stress tensor, and the two should not be 
confused. 

An outline of this work is as follows. Section 2 presents a simple discussion of cor- 
related fluctuations and mean flow, in particular emphasizing the circumstances in which 
fluctuations, not just mean flow values, must be retained in the governing equations. Sec- 
tion 3 focuses on turbulent energy transport. The key result is equation (15), which is the 
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fundamental thermodynamic relation for the fluctuations. Section 4 develops simple, one- 
dimensional (ID) solutions to the governing equations, contrasting them with solutions from 
the literature. Finally, §5 is a summary of the paper. 



2. Correlated Fluctuations and Mean Flow 



As stated in the Introduction, our presentation is based on weak turbulence theory, 
which bridges the divide between phenomenological and first-principle approaches to tur- 
bulent flow. Any flow quantity may be decomposed into the sum of a mean value plus a 
fluctuation with zero-mean. We use the notation {X) for the mean value of quantity X, and 
5X to denote its fluctuating component. For the flow velocity, weak turbulence means that 
the fluctuations are small compared to the isothermal sound speed, cs- For the pressure 
and density, the fluctuations are to be considered small compared with the mean value. The 
magnetic field B is assumed to be subthermal, i.e., the Alfven speed v a ~ Sv <C cs- In what 
follows, we shall not need to distinguish between the mean and fluctuating components of 
B. Mean velocity components may be either very large or very small compared with their 
RMS fluctuations. The rotational velocity will generally greatly exceed its azimuthal veloc- 
ity fluctuation, but the mean radial drift velocity, a consequence of dissipation, will be of 
second order (i.e., quadratic) in the fluctuation amplitude. 

This last point is key, and merits some discussion. Upon long term time-averaging, the 
governing equations of a turbulent flow have a leading asymptotic form which is either zeroth 
or second order in the amplitudes. It is important to distinguish these two cases, because 
by definition fluctuations play no role in the leading behavior of a zeroth order equation, 
whereas fluctuations must be self-consistently retained everywhere in a second order equation. 
Consider first the exact radial equation of motion expressed in the usual (R, 0, Z) cylindrical 
coordinates: 

d _\ vl idP d<5> 1 



m+v . V ) VR ^ = -- m - m - + - (JxB)R . (i) 

Our notation is standard: v is the velocity (subscripts denote components), p the mass 
density, P the gas pressure, $ the external gravitational potential, and J the current density. 
If the mean radial velocity is of second order and the magnetic field is weak, the time-steady 
leading order form of this equation is simply 

Rtf - 1 9{P) + ^ (2) 

where Q is the mean angular velocity and all other terms are leading zeroth order mean 
values. In accord with familiar treatments of accretion flows, fluctuations play no role here. 
(A similar argument can be made for the z equation of hydrostatic equilibrium.) Contrast 
this, however, with the azimuthal equation of motion: 

Here, there are no leading terms of zeroth order: after averaging, the leading behavior of 
the equation is second order in the amplitudes (Balbus & Papaloizou 1999). All terms in- 
volve either the second order radial drift velocity (vr), or explicit amplitude correlations of 
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the form (SXSY). These latter turbulent correlations must be retained, or the resulting 
equation is incorrect. Indeed, the correlated fluctuations in this equation are just the (non- 
advective) angular momentum flux, which is critical to the accretion process itself. Standard 
a treatments in effect retain these correlations by inventing an "anomalous viscosity." To 
the extent that the theory is independent of the functional form of the correlations, the fact 
that at least something is present in the equation where the turbulent stress should appear 
is sufficient for some purposes. It allows points of principle to be illustrated (e.g., the initial 
spreading and subsequent concentration of an evolving disk), as well as some robust results 
to be obtained (e.g., the luminosity-accretion rate relation [Pringle 1981].) 

Genuine difficulties with the standard treatment emerge when one attempts to go beyond 
rotational mechanics. Consider the entropy equation: 

t^i (it + v ■ v ) lnPp ~ 7 = Q+ ~ g ~' (4) 

where Q + and Q~ are respectively the volume specific entropy gains and losses. The tacit 
assumption that is almost always made is that the flow quantities appearing in equation (4) 
may be replaced by their mean values. This is incorrect. 

The problem is that the mean radial velocity is a second order quantity, and the entire 
entropy equation is therefore also of second order. The quantity Pvr is not, on average, 
simply the product of the mean pressure and radial drift velocity, but (in one dimension): 

(Pvr) = (pOv R ) = (9)(pv R ) + (p)(5v R 59) + higher order terms, (5) 

where 9 = c| is the temperature variable. In fact, the full entropy equation is even more 
complex, because there are correlations between 8v R and the fluctuation of the entropy 
gradient term itself that must be taken into account. There is no a priori reason that all of 
these correlation coefficients should vanish. Clearly, care must be taken when determining 
quantities such as the sign of the mean entropy gradient in a turbulent accretion flow. A 
more systematic approach is needed. 



3. Turbulent Transport in NRAFs 

3.1. Second Order Conservative Transport Equations 

To investigate the behavior of a model flow with a turbulent thermal energy flux in more 
detail, begin with the time-steady form of the angular momentum conservation equation, 



R 

V- ( pRv+v - ^ B <p B , 







(6) 



and the energy conservation equation: 
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We have neglected the contribution of the particle viscosity and (potentially more seriously) 
the thermal conductivity. Throughout the remainder of this paper, all equations are un- 
derstood to be azimuthally averaged. Unsubscripted boldface vectors and vector operators 
should be regarded as poloidal. 

We next expand all quantities into mean plus fluctuating components. The largest 
terms in the energy flux are second order in the fluctuating amplitudes. (Recall that we are 
assuming that the Alfven speed is of the same order as a kinetic velocity fluctuation.) The 
mass flux is 

(pv) = (p)(v) + (5 P 8v), (8) 

and satisfies 

V-(pw) = 0. (9) 
If we now define the stress vector W to be 

W=(5v*8v-?j£), (10) 

then the angular momentum equation becomes 

V- ((pv}R 2 n + (p)RW) = 0. (11) 
The energy equation, through second order in weak turbulence theory, reads: 



R 2 W , „ , 1(9) \ , , , 7 



(pv) [ — + $ + -i^Lj + (p)mW + -J—(p)(508v) 



-Q- (12) 



The first group of terms multiplying the mass flux (pv) corresponds to the Bernoulli con- 
stant in polytropic spherical flow. In general, of course, it need not be constant. The key 
component of the energy flux that distinguishes transport in a nonradiative flow from a thin 
Keplerian disk involves only one new correlation product: the thermal energy flux propor- 
tional to (S65v). This term is generally ignored in a models of nonradiative accretion flows. 
Instead, these models retain a form of the energy flux appropriate for a radiative thin disk, 
in which the dominant nonadvective energy flux comes from the rotational transport term 
W. In general, the ratio of the thermal energy flux to rotational stress will be, as noted, 
of order cs/RQ- While this justifies neglect of the thermal energy flux in thin disk models, 
it also shows that this correlation must be retained in any flow with comparable rotation 
and sound speed. Indeed, retention of the energy flux terms is essential to formulating a 
thermodynamically self-consistent model for the energetics of turbulent fluctuations. 



3.2. Fluctuation Thermodynamics 

The result of combining mass conservation (9) with radial hydrostatic equilibrium (2) in the 
energy conservation equation is 
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where 

S = \n[(P){p)-i]. (14) 

If we now combine this equation with angular momentum conservation (11), our result 
simplifies to 

v. (^(fc»>) + ^W-vs = -(«- + (is) 

Equation (15), which is the internal energy equations for the fluctuations, is a key the- 
oretical result of this paper. It is readily interpreted. The right hand side is the rate at 
which fluctuations exchange energy with their "surroundings." The first term (Q~) repre- 
sents bulk radiative losses, and the final term is the rate at which energy is supplied by the 
free energy reservoir of differential rotation. In a classical thin Keplerian disk, these terms 
together comprise the dominant energy balance, and the entire left side of the equation may 
be ignored. 

In a nonradiative flow, however, none of the free energy of differential rotation escapes, 
but is instead routed through two possible thermodynamic paths. The first possibility, 
embodied by the divergence term on the left side of (15), is direct mechanical transport of 
the thermal energy content of a fluctuation, including simple advection plus PdV work. The 
form of the flux is 

l@-(6e8v) = (p)Cp(5T6v), 
7-1 

where Cp is the mass specific heat capacity at constant pressure (cf. Schwarzschild 1958). In 
the case of an adiabatic sound wave, this expression is equivalent to (SPSv), precisely the 
formal rate of work done by the wave. The balance set by equating this term with the final 
term on the right side of equation (15) corresponds to wave action conservation for adiabatic 
disturbances (Lighthill 1978). 

The second term on the left is the average rate at which entropy is generated by the 
mean mass flow. It is the only other thermodynamic channel available to the free energy 
of differential rotation: if the energy is not carried off by adiabatic processes, then it must 
be dissipated. In turn, the dissipated energy must, in the absence of radiative losses, drive 
either an inflow or outflow. Dynamical considerations strongly favor an outflow, at least in 
the outer flow regions, leading to the scenario envisaged by Blandford & Begelman (1999). 
To the extent that bulk flow does not occur, a thermal flux must carry the energy away. This 
term does not appear in the Blandford & Begelman formulation, so this form of transport 
is precluded in their analysis. 

Finally, it should be noted that equation (15) differs somewhat from the internal equa- 
tion (4) used in Abramowicz, et al. (2002) and elsewhere by the same authors. This is 
because the fluctuation formalism is not used by Abramowicz, et al. (2002), so the internal 
energy input is regarded as a source of particle heating, and is therefore dissipative. By 
contrast, in this work we have granted the fluid the degrees of freedom associated with fluc- 
tuations, which allows the free energy of differential rotation to be extracted with or without 
irreversible dissipation. 
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3.3. Thermal Convection vs. Thermal Flux 

The discussion of the previous section invites comparison with the convection-dominated 
flows put forth Narayan, Igumenshchev, & Abramowicz (2000). While the dynamical de- 
scription of these models has been criticized, (Balbus & Hawley 2002; Narayan et al. 2002 
for a reply), one should not dismiss the possibility that a thermal energy flux might affect 
the flow structure under more general circumstances. The point goes beyond the presence 
or absence of thermal convection, and allows us to sharpen the very notion of what it means 
for a flow to be "convection-dominated." 

In an accretion flow whose primary source is the free energy of differential rotation, 
what does it mean to say that the turbulence is dominated by thermal convection? The 
unstable linear mode in question is always the same slow mode branch of the MHD dispersion 
relation regardless of the mixture of adverse shear and thermal gradients that may be present. 
(Recent claims to the contrary in the literature reflect confusion on this point.) If the primary 
instability is the MRI, which in its nonlinear resolution allows accretion to proceed, and any 
adverse entropy gradients are present only as a secondary consequence of accretion, in what 
sense can the flow ever be dominated by convection? 

As originally formulated (Narayan, Igumenshchev, & Abramowicz 2000; Abramowicz, 
et al. 2002), the role of convection was truly dominant: it was argued that it could cause 
the stress tensor to vanish and that it could eliminate turbulent dissipation throughout the 
body of the accretion flow. But this leads to serious difficulties (Balbus & Hawley 2002). 
The vanishing of the radial component of W (more precisely a density-weighted average 
thereof, denoted T R ^) leaves the fluctuations without a local source of free energy. Source- 
free, dissipation-free local turbulence maintained by large dissipative entropy gradients is 
a thermodynamically dubious proposition. It is therefore not surprising that global, three- 
dimensional MHD accretion simulations (Hawley, Balbus, & Stone 2001; Hawley & Balbus 
2002; Igumenshchev, Narayan, & Abramowicz 2003) all report robust positive values of Tr^. 

More recently, claims that a simulated accretion flow is dominated by convection have 
shifted to morphological grounds (Igumenshchev, Narayan, & Abramowicz 2003), with 
only the visual appearance of the flow used as a basis for an identification. This is not 
the diagnostic of choice. Whether or not thermal convection is of significance in a high 
temperature NRAF simulation is in fact a well-posed quantitative problem. One must first 
begin by treating energy and angular momentum transport on the same footing by elevating 
the energy correlation flux (595v) to the same status as T R ^. Second, controlled studies are 
essential. The visual appearance of turbulence is not a very reliable guide to the presence or 
absence of thermal convection, especially if the root cause of the turbulence is a magneto- 
rotational process. Simulations must be run with and without a dissipative heat source, 
and both turbulent correlation coefficients examined in detail for each case. This has in fact 
already been done for Tr^: Stone & Pringle (2001) and Hawley & Balbus (2002) both found 
very little difference in the behavior of the stress tensor whether or not resistive heating was 
present in their simulations. The turbulent energy flux has yet to undergo a similar level 
of scrutiny, however. The question is, how much does this thermal flux correlation change 
relative to its value in the absence of a dissipative heat source? What matters ultimately 
is whether the thermal energy flux is an efficient drain for the free energy of differential 
rotation in any nonradiative flow, convective or otherwise. 
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4. One-Dimensional Solutions 

4.1. Formulation of the Problem 

To better understand the physical consequences of a thermal energy flux in an NRAF, we 
apply the formal results of §3 to some explicit, simple accretion solutions. To do this, we 
need explicit functional forms for W and (59dv). Since our goal is one of illustrating a point 
of principle, we will use the dimensional analysis associated with a formalism for each of 
these fluxes. That is, 

{pSvuSv* - = ass(P), (p5v R 59) = a T (P)(c s ) (16) 

where ass (the Shakura-Sunyaev a parameter) and «t are dimensionless constants. Note 
that ass > 0, but we leave the sign of «y unconstrained for the moment. The dynamical 
stress must be positive in order to extract energy from the differential rotation, whereas the 
the thermal flux is a secondary consequence of the turbulence, and could in principle have 
either sign, depending on what mode of transport dominates. This is particularly important 
when one uses this flux as a numerical diagnostic, because the restrictions sometimes used 
(e.g. two dimensions, adiabatic gas law) can bias the direction of the energy flux. 

Next, we take the rather drastic but standard step of height-integrating the equations 
of motion, and assume that what remains is a well-defined, ID (on average) radial flow 
(Narayan, & Yi 1994; Blandford & Begelman 1999). It is useful to have such a formal 
solution to compare with others in the literature. Numerical NRAF simulations do show a 
well-defined thick wedge forming along the midplane with a relatively autonomous interior, 
surrounded by a low density coronal envelope (Hawley & Balbus 2002). There is some sense 
in regarding much of the wedge structure as ID, although there is also a significant amount 
of high latitude backflow that is lost in this approximation. Perhaps the best reason to study 
these ID flows is to illustrate a point of principle when they fail to work self-consistently, in 
contrast, say, to steady viscous disk models. 

Our ID formulation is similar to the "idealized ADAF" analyzed by Blandford & Begel- 
man (1999), with an added thermal energy flux. The equations of motion, which can now 
be expressed in terms of mean value quantities, are 



trifl . „ Cj 



2l +a - p = i? (18 > 

m fR 2 n 2 GM 1 V\ „,„ „ R't V 3 ' 2 „ 

These are, respectively, the radial equation of motion, and the angular momentum and 
energy conservation equations, m is the conserved mass accretion rate, Cj is the conserved 
rate of angular momentum transport, and Ce is the conserved energy transport rate. S 
is the height integrated density (i.e., the column density), and V is the height integrated 
pressure. We have assumed a central black hole mass M, but as is customary, have used a 
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gravitational potential depending only upon R. An alternative useful form for equation (19) 
is obtained by eliminating the ass term via equation (18): 

m fR 2 n 2 GM 7 P\ R 1 V 3 ' 2 _ ^ 

^{— + nr-—z) + — aT ^ = CE - nCj (20) 



4.2. Static Envelope 

We begin with a search for solutions with fa — 0. Combining equations (17) with (18) gives 

Cj 



a ss R 2 ' 

and 



(P) oc R-', (21) 



n2 GM 3 Cj 

U ~ R 3 O^Ei^ 1 j 

The radial power law dependence of the pressure suggests we look for similar behavior in fl 
and E. This is possible only if Ce = 0, so that the net energy flux vanishes. Then equation 
(20) leads directly to 

n = -—^iw^- (23) 

Evidently, it is necessary that ax < 0, i.e. the heat transport must be inward to ensure there 
is no energy deposition. Eliminating fl between the last two equations leads to an equation 
for E: 

and Q follows from (23), 

n» = (25) 

1 + ^ R 3 K J 

This is a sub-Keplerian profile and is associated with thick disk structure (Frank, King, & 
Raine 2002). The temperature is given by 

W = |=™ (1 + { )-. (26) 

Equations (24-26) comprise our ID static solution. Notice that it satisfies the virial equilib- 
rium condition 

2T + $ + 3(#) = (27) 

where T = R 2 Vt 2 /2 is the rotational kinetic energy, and that it differs from the static convec- 
tion dominated solutions of Narayan, Igumenshchev, & Abramowicz (2000) and Abramow- 
icz, et al. (2002) in requiring an inward, rather than an outward, thermal flux. The cause 
of this difference is that our solution has a positive value of the R<fi stress tensor component, 
as is required by energy conservation. 
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4.3. Accreting Envelope 

Equations (17-19) allow for power law solutions even when rh does not vanish, a property 
shared with NRAFs in their standard formulation (Narayan, & Yi 1994; Blandford & Begel- 
man 1999). We require the vanishing of both the energy as well as the angular momentum 
flux, Ce = Cj = 0. Since Wj^ is no longer directly proportional to Cj, as in the static solu- 
tion of the previous section, there is no difficulty with simultaneously demanding a positive 
stress and a vanishing net angular momentum flux. This condition is satisfied in classical 
Keplerian disks at radii large compared with the inner edge. 

With 

VocR- 3 / 2 , Soci?" 1 / 2 , fioc/T 3 / 2 , (28) 

equation (17) becomes 

V 2 (GM » 9 \ . . 

- R 2 Q 2 (29) 



E 5 \ R 

and after some manipulations, the energy equation (19) may be written 
Defining 

T B = 1 ~T> C = 7 ^UJ ' X = GM' (31) 

we find that x must satisfy 

x\C 2 + A 2 /4)-x 2 (C 2 + AB) + B 2 = 0. (32) 

The choice of sign for this quadratic (in x 2 ) equation must be chosen to be consistent with 
equation (30), a spurious root having been introduced in the course of squaring a radical. 

It is the appearance of C that distinguishes our solution from those of earlier NRAF 
studies. If C > 0, corresponding to outward transport, then it must be large compared 
with A to change the C = solution x 2 = 2B/A qualitatively. (The novel B — solution 
of equation [30] is spurious when C > 0.) On the other hand, if C < 0, corresponding to 
inward transport, a qualitative change does occur: whereas a standard ID NRAF has no 
solution for the important case 7 = 5/3 (i.e. B = 0), the presence of an inward energy flux 
leads to a new sub-Keplerian solution 



x 



2 



'l + A 2 /AC 2 )-\ 9 = (2/h)(GM/R)(l + AC 2 /A 2 )- 1 . (33) 



An inward thermal flux seems anomalous. Considering also its novel consequences for both 
static and accreting flows, we are thus led to consider the physical processes leading to either 
inward or outward transport. 
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4.4. Discussion 

It is perhaps surprising that, in these ID models, inward thermal transport is required to 
establish either a static profile, or a 7 = 5/3 accretion profile. But given a power law leading 
asymptotic order behavior for the thermal and rotational energy fluxes at large R, there is 
little choice. Energy must be strictly conserved. Therefore, the two fluxes must either cancel 
one other, or their sum must be conserved. Because it forces a unique power-law behavior, 
the latter possibility overconstrains the problem. Hence the rotational energy flux, which 
must always be outward, is canceled by an equal and opposite thermal flux. 

Second, equation (15), together with static or isentropic flow, implies a very frugal use 
of free energy: essentially all of it must go into an inward thermal flux. What might be the 
physical basis for such behavior? 

The MRI, acting in the presence of a Schwarzschild-stable entropy gradient, will tend to 
drive a thermal energy flux down the entropy gradient, which in this case is indeed inwards 
(Balbus 2000). For the static solution, equations (21) and (24) give 

S ~ \nR 2 ^ 3 

for the static solution, which is Schwarzschild-stable for 7 > 1.5. On the other hand, for the 
7 = 5/3 accreting solution, S is constant with R. 

But there are deeper problems with the mechanism of turbulent mixing, even for the 7 > 
1.5 static solution. The inward energy transport is produced by an intrinsically dissipative 
process, and the entropy of a fluid element would continuously increase. Without radiative 
losses, a static solution would not persist under these circumstances. This is reflected in 
equation (15), where it may be seen that the dissipative term is coupled both to mass flow 
as well as to a finite entropy gradient. 

Might it be possible to continuously extract the free energy of differential rotation and 
have it go directly into nondissipative energy transport? Certainly. This is what waves do in 
the process of conserving their action, and non-dissipative wave propagation does not cause 
systematic mass flow in the background medium. But in our problem, any waves that are 
present will be coupled to dissipative turbulence. 

This raises another issue. In general, the linear behavior of a rotating magnetized 
fluid is marked by two incompressible modes. At a fixed wavelength and arbitrarily weak 
field strengths, one branch is a destabilized slow mode (the MRI); the other is generally a 
stable inertial wave (assuming that entropy gradient is dominated by its vertical component), 
which becomes Alfvenic as the field strengthens. Nonlinear coupling between the unstable 
and stable modes is almost certain to be present in a turbulent fluid. 

Hydrodynamical inertial waves in disks are characterized by frequencies u below the 
epicyclic value (Vishniac & Diamond 1989; Balbus 2003), which in a Keplerian disk is 
simply Q. Hence axisymmetric waves propagate inwards, because any initially outward 
propagating wave at fixed cu will eventually encounter a turning point at which u exceeds k, 
and the wave must reflect. Two-dimensional numerical simulations of the MRI would then 
be characterized by ass > (outward transport by rotational stress) and, to the extent that 
inertial forces exceed those of magnetic tension, c*t < 0. The MRI dominates the angular 
momentum transport (axisymmetric inertial waves carry no angular momentum of course), 
but perhaps not the energy transport. A preliminary numerical investigation of a global, 
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two-dimensional, self-gravitating, adiabatic, magnetized torus does indeed show comparable 
values of ass and \cxt\, with cut < (Fromang, Devilliers, & Balbus 2003). 

In three dimensions, the situation is significantly more complex, because inertial waves 
are not free to propagate inwards without reflection. Rather than uj < k, the condition for 
propagation is 

\uj — mQ\ < k, 

where m is the azimuthal wave number. The radius at which equality in the above relation 
is obtained is known as the Lindblad resonance, and normally there are two such locations: 
outer and inner Lindblad resonances, where u — mQ is respectively positive and negative. 
The existence of an inner Lindblad resonance is a consequence of nonaxisymmetry, and the 
resulting trapping of any inertial waves undermines the arguments used in the axisymmetric 
case. One may always appeal to the magnetic field to get passed the Lindblad barriers, but 
there is no guarantee that waves liberated in this way will be heading preferentially inwards. 

The case for inward thermal energy transport is even more problematic in real black hole 
systems. Coulomb conduction, heretofore ignored, is likely to be an important process in the 
dilute plasmas that characterize nonradiative black hole accretion. The relevant instability 
criterion for a plasma whose thermal conductivity is dominated by a magnetized Coulomb 
conductivity is that the temperature, not the entropy, increase outward (Balbus 2001). 
The combined magnetorotational/magnetothermal instability would result in outward, not 
inward, thermal transport. A static solution under these conditions is not possible. 

Both acoustic waves and turbulence associated with an unstable thermal gradient will 
tend to move thermal energy outward with high efficiency, since in both cases the correlation 
between the velocity and temperature fluctuation is strong. In the case of acoustic waves, a 
net outward flux is expected since initially inward propagating waves will tend to refract and 
head outwards. Outwardly propagating waves, on the other hand, will become increasingly 
dominated by their radial wavenumber, and eventually dissipate or shock. This would seem 
more likely to result in a wind from the outer regions of the flow, rather than a static or 
accreting halo. 

To conclude, in a height-integrated ID analytic formalism, inward transport of thermal 
energy is needed for a static solution or for 7 = 5/3 rotating accretion to exist. Conditions fa- 
vorable to inward transport include axisymmetric and adiabatic flow, restrictive assumptions 
likely to be more common in numerical simulations than in nature. Three-dimensional flow, 
and the mixing of magnetothermal with magnetorotational instabilities in the presence of a 
Coulomb conductivity, all tend to produce an outward energy flux. In the presence of a pos- 
itive rotational stress, outward thermal transport does not stifle accretion. It does, however, 
lead to accretion flows whose salient properties cannot be captured by ID height-integrated 
modeling. 



5. Summary 

The key point of this paper is that thermal fluctuations must be taken into account in the 
formulation of governing equations of nonradiative accretion flows, and in particular that 
the entropy equation cannot be analyzed directly by assuming all flow quantities take on 
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their mean values. The energy flux 
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is an important component in any nonradiative accretion process. Energy conservation, the 
central assumption of nonradiative flows, is not satisfied unless this term is present, and 
such basic processes as wave action conservation are lost. It is important here to note that 
the presence of this correlation flux is more fundamental than the mere inclusion of a some 
additional transport process, such as thermal conduction. The issue is one of thermodynamic 
self-consistency. The energy flux (34) must accompany the advected flux term 



for the same reasons, and at the same level of formalism, that the stress tensor W accom- 
panies the advected angular momentum flux. The results are summarized quantitatively in 
equation (15), which expresses energy conservation for either wavelike or turbulent fluctua- 
tions. 

With one exception, no other accretion model familiar to the author includes a non- 
advective thermal energy flux. As noted in §3.1, convection-dominated models (Narayan, 
Igumenshchev, & Abramowicz 2000; Quataert & Gruzinov 2000; Abramowicz, et al. 2002) 
assume that the energetics are dominated by an outward energy flux generated by a convec- 
tively unstable accretion process. The central assumption of these models is that convection 
inhibits angular momentum transport and causes the turbulent stress to vanish. We have 
argued, on the other hand, that the stress tensor Wj^ must be positive in order to extract 
free energy from differential rotation. One consequence of this is that any viable, height- 
integrated static solution must have an inward flux of thermal energy, rather than an outward 
flux, as expected in a convectively unstable flow. (The same inward flux requirement holds 
for the existence of steady, height-integrated, 7 = 5/3 rotating accretion flows.) It there- 
fore seems unlikely that thermal convection, even if it were present, would stifle accretion. 
Instead, it would simply add to the outward rotational energy transport while negligibly 
affecting the angular momentum transport (Stone & Pringle 2001; Hawley & Balbus 2002). 

Establishing how the thermal flux emerges in a mean flow formalism for weak turbulence 
is a theoretical result that we believe will foster a deeper understanding of nonradiative 
accretion, if for no other reason then for avoiding the pitfall of confusing flow quantities 
with their mean values. But the immediate utility of the present work is likely to be as a 
numerical diagnostic. This is particularly true for those cases in which thermal convection is 
claimed to be present, to the Reynolds/Maxwell stress tensor, the energy correlation {6v 59) 
has not been well-studied. It is particularly important to know how much of the free energy 
is ultimately "radiated" by this under-explored form of mechanical luminosity. 
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